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摘 要: 基于 1954 一 2018 年 皇甫 川 流 域 长 时 间 序 列 降水 量 .径流 量 和 输 沙 量 监测 数据 ,运用 线性 趋势 法 Mann- 
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Kendall 检验 .累积 距 平 法 等 方法 研究 了 皇甫 川 流 域 水 沙 周期 演变 特征 ,应 用 水 文 分 析 法 定量 评估 了 降水 变化 和 人 
类 活动 对 皇甫 川 流 域 径流 量 和 输 沙 量变 化 的 贡献 率 。 研 究 结果 表 明 : 旦 甫 川 流 域 降水 量 随 时 间 变 化 旦 不 明显 减 小 
趋势 ,而 径流 量 和 输 沙 量 随时 间 变 化 却 呈 显著 减 小 趋势 , 且 径 流量 减少 趋势 大 于 输 沙 量 减 少 趋势 。 径 流量 和 输 沙 


量 均 从 1986 年 发 生 巾 丰 到 枯 的 突变 。 以 1954—1986 年 为 基准 期 ,1987 一 2018 年 人 类 活动 对 皇甫 川 流 域 径流 量 和 


输 沙 量 的 影响 分 别 为 83.04% 和 75.24% ,降水 变化 对 皇甫 川 流 域 径流 量 和 输 沙 量 影响 分 别 为 16.96% 和 24.76%。 突 


变 年 前 后 输 沙 量 出 现 明显 的 下 降 , 产 流 能 力 对 输 沙 量变 化 的 贡献 率 超过 了 85%。 降 水 变化 和 人 类 活动 是 造成 皇甫 
川 流 域 水 沙 减 少 的 影响 因素 ,其 中 流域 内 逐年 增强 的 人 类 活动 是 促使 皇甫 川 流 域 径流 和 输 沙 减少 的 主要 因素 ,而 


气候 变化 引起 的 降水 减少 是 次 要 因素 。 
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河流 沉积 物 经 河流 从 大 陆 到 海洋 的 输 移 是 调 
节 河 岸 稳定 .土壤 形成 .元 素 生物 地 球 化 学 循环 号 、 
地 壳 演 化 等 与 地 球 有 关 的 重要 过 程 之 一 ” ,然而 , 随 
着 全 球 变化 加 剧 ,可 利用 水 资源 正在 枯竭 , 输 向 海洋 
的 泥 沙 量 正在 减少 ? ,给 流域 管理 带 来 许多 挑战 后 。 
黄河 是 著名 富 含 泥 沙 的 河流 ,其 流 经 平原 地 区 由 于 
河道 拓宽 .坡度 减 小 ,导致 黄河 中 上 游 大 量 输 沙 在 
下 游 沉积 外 ,再 加 上 气候 变化 和 人 类 密集 活动 E 
地 利用 变化 .河流 管理 .兴修 梯田 、 筑 坝 、 退 耕 还 林 
等 ) ,导致 径流 量 和 泥 沙 量变 化 ,而 这 种 变化 关系 到 
整个 黄河 流域 生态 稳定 "。2000 年 以 来 ,黄河 流域 
输 沙 量 出 现 显著 减少 趋势 ,而 这 种 变化 受到 国内 
外 学 者 的 关注 ,认为 黄土 高 原 地 区 “退耕 还 林 ( 草 )” 
工程 实施 ,逐年 增加 了 植被 覆盖 度 UE ROMPE E 
表 径 流 ,减少 了 土壤 侵蚀 ,调控 了 河道 水 沙 , 使 得 黄 
土 高 原 地 区 生态 环境 大 大 改善 ,水 沙 变化 情势 发 生 
了 新 的 变化 ”。 

皇甫 川 是 黄河 的 一 级 支流 ,在 20 世纪 50 年 代 
每 年 向 黄河 输送 泥 沙 0.64x10: (5 ,严重 的 水 土 流失 
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从 20 世 纪 60 年 代 以 来 ,皇甫 川 流 域 实施 了 一 系列 
水 土 保持 措施 ,包括 梯田 、 大 坝 \ 水 库 建 设 .退耕 还 
林 、 植 被 恢复 等 ,造成 流域 径流 量 和 输 沙 量 开始 减 
少 ,改变 了 黄河 下 游 水 沙 情 势 ”…。 新 时 期 皇甫 川 流 
域 水 沙 将 如 何 变化 ?引起 流域 水 沙 关 系 变 化 的 内 
在 驱动 因素 又 是 什么 ”这些 问 题 将 会 对 皇甫 川 流 
域 生态 可 持续 发 展 和 水 沙 调 控 宏观 决策 产生 重要 
影响 。 基 于 此 ,人 研究 采用 长 时 间 序 列 水 沙 资料 对 流 
域 径 流 和 输 沙 演变 规律 进行 分 析 ,探究 降水 变化 和 
人 类 活动 对 皇甫 川 流 域 径 流量 和 输 沙 量变 化 的 页 
献 率 ,以 期 为 皇甫 川 流 域 水 土 保持 建设 ,生态 环境 
恢复 ,建立 流域 水 沙 科学 的 调控 体系 提供 理论 依据 
和 决策 参考 。 


1 研究 区 概况 与 方法 


1.1 研究 区 概况 
£r JI] RUE ( 39912' -39954'N, 110918' -112912/E) 
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位 于 黄河 中 游 河 口 镇 与 龙门 区 间 右 岸上 段 , 属 于 黄 
河 的 一 级 支流 , 源 于 内 蒙古 东部 鄂尔多斯 准 格 尔 旗 
点 上 畔 沟 , 流 经 黄土 高 原 过 渡 带 和 内 蒙古 草原 ,在 陕 
西 省 府 谷 县 的 巴 免 坪 汇 人 黄河 ,流域 总 面积 3246 
km ,有 一 级 支流 49 条 ,其 主要 干流 纳 林 川 全 长 137 
km ,流域 面积 1351 km。 皇 甫 川 流 域 气候 类 型 为 温 
带 大 陆 性 季风 气候 ,年 平均 气温 9.1 ,降水 主要 集 
中 在 6 一 9 月 ,降水 量 占 全 年 的 80% 以 上 ,多 年 平均 
降水 量 为 350~450 mm ,气温 和 降水 在 空间 分 布 上 具 
有 从 东南 向 西北 地 区 逐渐 减少 的 变化 特征 。 流 域 
沟 蜜 纵横 ,水 土 流 失 严 重 , 土壤 侵 刨 以 水 力 侵蚀 为 
主 ,并 复合 有 风力 .重力 侵蚀 。 流 域 地 貌 具 有 梁 往 
FZI 沟 穴 密布 .地形 破碎 的 特点 ,黄土 厚度 由 西北 
向 东南 递增 ,土壤 结构 疏松 , 极 易 风化 , 遇 高 强度 暴 
雨 极 易 产 生 侵蚀 ,其 现 有 的 植被 类 型 主要 为 草地 、 
农田 和 灌木 林 和 乔木 林 , 过 去 分 布 的 天 然 林 和 灌木 
草原 大 部 分 已 经 被 人 工 植被 所 代替 。 
1.2 数据 来 源 

考虑 到 资料 连续 性 和 完整 性 ,根据 流域 内 水 文 
站 点 分 布 ,研究 基于 1954 一 2018 年 皇甫 川 流 域 径流 
泥 沙 数 据 来 分 析 水 沙 周 期 变化 特征 ,其 中 1954 一 
2009 年 径流 和 泥 沙 数据 来 自 于 黄河 流域 水 文 年 鉴 ， 
2010—2018 年 来 自 于 黄河 泥 沙 公报 。 皇 甫 川 境内 
有 雨量 站 7 个 ,但 是 部 分 雨量 站 数据 有 缺失 ,采用 将 
7 个 雨量 站 的 数据 进行 反 距离 加 权 后 形成 的 1954 一 
2018 年 皇甫 川 流 域 降水 量 数据 ,用 于 讨论 气候 变化 


-一 河流 
CO 黄河 流域 


图 1 皇甫 川 流 域 地 理 位 置 示意 图 
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对 皇甫 川 水 沙 变 化 的 影响 。 根 据 一 些 相关 研究 成 
RE ,黄土 高 原 气 温 与 径流 量 和 泥 沙 量 没 有 显著 的 
相关 性 ,所 以 在 讨论 气候 变化 对 流域 水 沙 影响 因素 
中 只 考虑 了 降水 因素 。 人 类 活动 影响 径流 量 和 输 
沙 量 有 多 种 方式 ,但 是 各 项 活动 很 难 根 据 以 年 为 单 
位 进行 分 项 统计 ,人 类 活动 对 水 沙 变化 影响 数值 很 
难 收集 ,加 上 不 同年 代数 据 统 计 标 准 不 一 致 ,造成 
一 些 统计 资料 失实 ,很 难 分 项 统计 各 项 人 类 活动 措 
施 对 水 沙 变化 影响 ,因此 ,研究 中 仅 讨 论 了 人 类 活 
动 对 皇甫 川 流 域 水 沙 变化 的 贡献 率 ””。 

1.3 研究 方法 

1.3.1 径流 量 和 输 沙 量 阶段 性 变化 ”通过 简单 的 线 
性 回归 分 析 , 运 用 累积 距 平 法 系统 分 析 皇 青川 径流 
量 和 输 沙 量 长 时 间 序 列 阶 段 性 差异 特征 。 对 于 特 
定时 间 序 列 庆 (YX, Xo, Xs, X), KEX n, 累积 距 平 
EAn h PST ARE: 


X," Y(X,-X) 1-12, n (1) 
i=1 


式 中 : 怠 为 时 间 序 列 观测 值 ; X 是 为 平均 值 ;XxX 代表 
第 i 年 累积 值 。 累 积 距 平 值 主要 检测 观测 值 可 能 存 
在 的 突变 点 ,以 及 观测 值 随时 间 序 列 变化 的 趋势 。 
当 观 测 值 为 正 值 时 ,说 明 突 变 点 的 值 高 于 平均 值 ， 
相反 ,为 负 值 则 说 明 突 变 点 的 值 低 于 平均 值 。 

1.8.0 输 沙 归 因 诊断 分 析 ” 输 沙 归 因 诊断 分 析 是 用 
于 判断 降水 、 产 流 能 力 和 产 沙 能 力 在 皇甫 川 输 沙 量 
比例 变化 率 上 的 贡献 ,从 而 判断 出 各 因子 对 输 沙 量 
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Fig. 1 Schematic diagram of the geographical location of the Huangfuchuan basin 
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FNS R=0.01, P>0.1 
s-P( E S)ePxess (2) 600 
PN E 500 
式 中 :P 为 皇甫 川 流 域 平均 降水 量 ;f 为 产 流 能 力 ( 径 400 
流量 与 降水 量 的 比值 ) ;* 为 产 沙 能 力 ( 输 沙 量 与 径 8 300 


流量 的 比值 );F 为 径流 量 ;5 为 输 沙 量 。 
该 方法 定义 X() 的 比例 变化 率 为 :r(X)=X'dX/di， 

则 输 沙 量 $ 的 比例 变化 率 为 : 
r(S) e r(P) v r(f) +r 


2 结果 与 分 析 


2.1 降水 量 年 际 变化 特征 

从 表 1 和 图 2 可 以 看 出 ,1954 一 2018 年 皇甫 川 
流域 年 均 降 水 量 为 419.78 mm ,变异 系数 为 0.30, 属 
于 中 等 尺度 变异 。 最 大 降水 量 为 715 mm, 出 现在 
1967 年 ,最 小 值 为 211.3 mm, 出 现在 1965 年 。20 世 
纪 皇 甫 川 流 域 各 代 际 降水 量 呈 逐步 减少 的 趋势 , 进 
和 21 世纪 降水 量 才 呈现 出 逐渐 增 大 的 趋势 。 降 水 
量变 异 系数 变化 范围 为 0.25~0.37, 其 平均 值 为 
0.30 ,属于 中 等 尺度 变异 。 降 水 量 极 值 比 介 于 2.31~ 
3.39 ,其 均值 为 2.73。 降 水 量 减少 速度 为 0.7lmm'.ar: 
(P»0.1) ,没有 达到 显著 性 水 平 ,表明 1954—2018 
年 皇甫 川 流 域 降水 量 呈 不 显著 减 小 趋势 。 


(3) 


R1 皇甫 川 流 域 降水 量 年 际 变 化 特征 
Tab.1 Annual variation characteristics of precipitation 


in Huangfuchuan basin 


变异 系数 极 值 比 
时 } 平均 值 / 2 

期 平均 值 /mm (C) (EVE) 
1954 一 1969 年 470.68 0.37 3.39 
1970 一 1986 年 398.72 0.25 2.47 
1987—2001 年 373.03 0.26 2.31 
2002 一 2018 年 444.68 0.29 2.75 
1954 一 2018 年 419.78 0.30 2.73 


表 2 皇甫 川 流 域 径流 
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图 2 皇甫 川 流 域 降水 量 年 际 变 化 趋势 
Fig.2 Linear trend of the Huangfuchuan basin for 


precipitation 


2.2 径流 量 和 输 沙 量 年 际 变 化 特征 
2.2.1 径流 量 年 际 变化 特征 ”从 表 2 和 图 3a 可 以 看 
出 ,1954 一 2018 年 皇甫 川 年 均 径 流量 为 1.23x10 
m’ ,变异 系数 为 0.74, 属 于 中 等 尺度 变异 ,径流 量变 
异 系数 变动 范围 为 0.56~1.01。 径 流量 极 值 比 介 于 
6.76~21 ,其 均值 为 14.51。 径 流量 变化 速度 为 -3.62x 
105m .ai(P<0.05) ,达到 了 显著 水 平 ,表明 1954— 
2018 年 皇甫 川 流 域 径流 量 呈 显著 减 小 趋势 。 
2.2.2 给 沙 量 年 际 变化 特征 ”从 表 2 和 图 3b 可 以 看 
出 ,1954 一 2018 年 皇甫 川 流 域 年 均 输 沙 量 为 3.7x 
10't, 变 异 系 数 为 0.97, 属 于 中 等 尺度 变异 。 输 沙 量 
变异 系数 介 于 0.63~1.39, 其 平均 值 为 0.97。 输 沙 量 
极 值 比 介 于 14.13~44.04, 其 平均 值 为 33.05。 输 沙 
量变 化 速度 为 -1.17x10; t-a” (P«0.05) ,达到 了 显 
著 水 平 ,表明 1954—2018 年 皇甫 川 流 域 输 沙 量 呈 显 
著 减 小 趋势 。 

径流 量 和 输 沙 量 趋势 检验 Z 值 均 为 负 , 通 过 了 
显著 性 为 99% 的 检验 ,表明 径流 量 和 输 沙 量 减少 趋 
势 明显 ,径流 量 趋 势 检 验 Z 值 为 -6.03 ,和 输 沙 量 趋势 
检验 Z 值 为 -5.76, 径 流量 减少 趋势 大 于 输 沙 量 减少 
趋势 。 径 流量 Sens 斜率 值 为 -0.036 , 输 沙 量 Sen's F} 


量 和 输 沙 量 年 际 变化 特征 


Tab. 2 Changes of average annual runoff and sediment load in the Huangfuchuan basin 


时 其 径流 量 输 沙 量 
平均 值 /10: m? 变异 系数 (C) ” 极 值 比 (EVR) 平均 值 /10*t 变异 系数 (C) ” 极 值 比 (EVR) 
1954 一 1969 年 2.07 0.64 12.36 0.61 0.80 32.75 
1970 一 1986 年 1.58 0.56 6.76 0.53 0.63 14.13 
1987—2001 年 0.92 0.75 17.92 0.29 1.07 44.04 
2002—2018 年 0.34 1.01 21 0.06 1.39 41.48 
1954—2018 年 1.23 0.74 14.51 0.37 0.97 33.05 
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(a) 径流 量 


y= -0.0362x*73.084 
R?-0.38 
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图 3 皇甫 川 流 域 径流 量 和 输 沙 量变 化 趋势 


Fig.3 Linear trend of the Huangfuchuan basin for runoff and sediment load 


率 值 为 -0.017, 径 流量 和 输 沙 量 Sen’s 斜率 值 也 都 为 
负 值 ,表明 皇甫 川 流 域 径流 量 和 输 沙 量 均 呈 下 降 趋 
势 , 且 径流 量 减少 趋势 大 于 输 沙 量 减少 趋势 。 
Mann-Kendall 趋势 检验 与 Sen's 斜率 值 的 研究 结 
一 致 。 
2.3 径流 量 和 输 沙 量 阶段 性 变化 特征 

从 图 4 可 以 看 出 ,皇甫 川 流 域 径 流量 和 输 沙 量 
累积 距 平 曲线 具有 明显 阶段 性 特征 。 皇 甫 川 径 流 
量 和 输 沙 量 累积 距 平 值 都 在 1986 年 达到 最 大 值 , 表 
明 径流 量 和 输 沙 量 都 在 1986 年 发 生 突变 。 根 据 突 


累积 距 平 输 沙 量 /10st 
oO — wt v0 5 Uu Oo a oo 


变 年 份 可 以 将 皇甫 川 流 域 径 流量 和 输 沙 量 分 为 
1954—1986 年 和 1987 一 2018 年 两 个 阶段 。1954 一 
1986 年 年 均 径 流量 为 1.84x10 mi ,1987 一 2018 年 年 
均 径 流量 为 0.66x10 mi ,突变 后 年 均 径 流量 比 突变 
前 减少 188%。1954 一 1986 年 年 均 输 沙 量 为 1.833x 
10 t, 1987—2018 4E. 4F H5] fig VP 5t 7 5.83 x 10 t, 突变 
后 多 年 年 均 输 沙 量 比 突变 前 减少 了 214%。 

根据 UF、UB 曲线 交点 位 置 ,可 以 确定 皇甫 川 径 
流量 和 输 沙 量 突变 都 是 从 1986 年 开始 (图 5a, 图 
5b) ,这 与 累积 距 平 检 验 结果 一 致 。 


累积 距 平 径流 量 /10s m? 


— 


年 份 


id 1 1 1 1 1 1 1 IN 0 
954 1961 1968 1975 1982 1989 1996 2003 2010 2017 


图 4 皇甫 川 流 域 径流 量 和 输 沙 量 的 累积 距 平 检验 


Fig.4 Testofcumulative anomaly of runoff and sediment transport in Huangfuchuan basin 
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图 5 皇甫 川 流 域 径流 量 和 输 沙 量 Mann-Kendall 突变 检验 


Fig. 5 Mann-kendall test of the Huangfuchuan basin runoff and sediment transport 
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2.4. 降水 变化 和 人 类 活动 对 皇甫 川 流 域 水 沙 变化 
的 影响 

以 皇甫 川 流 域 径流 量 和 输 沙 量 突变 点 以 1986 
年 为 分 界 点 ,将 皇甫 川 流 域 的 降水 量 、 径 流量 以 及 
输 沙 量 划 分 为 两 个 时 期 1954 一 1986 年 和 1987 一 
2018 年 。 

从 图 6 看 出 ,降水 -径流 与 降水 - 输 沙 双 累 积 曲 
线 都 出 现 了 拐点 ,突变 年 之 后 曲线 斜率 小 于 突变 年 
之 前 的 曲线 斜率 。 突 变 之 前 累积 径流 量变 化 幅度 
大 于 累积 输 沙 量变 化 幅度 ,累积 径流 量 曲线 斜率 大 
于 累积 输 沙 量 曲线 斜率 ,这 主要 由 于 降水 量 增加 ， 
使 河道 径流 量 增加 ,导致 河道 输 水 能 力 大 大 提高 ， 
并 且 这 种 能 力 大 于 河道 输 沙 能 力 ;突变 年 之 后 累积 
径流 量变 化 幅度 小 于 累积 输 沙 量变 化 幅度 , 累积 输 
沙 量 曲线 斜率 大 于 累积 径流 量 曲线 斜率 ,这 主要 由 
于 径流 量 减少 导致 输 沙 能 力 降低 ,大 量 泥 沙 汶 积 在 
河道 里 。 


将 突变 年 后 降水 资料 带 入 双 累 积 曲 线 建立 的 
回归 方程 ,得 到 降水 变化 和 人 类 活动 对 水 沙 变 化 影 
响 的 相关 关系 方程 式 , 其 相关 系数 达到 了 0.9 显著 
性 水 平 以 上 , 显著 性 水 平 忆 <0.05, 达 到 了 显著 水 
平 。 以 1954 一 1986 年 为 基准 期 ,1987 一 2018 年 人 类 
活动 对 皇甫 川 流 域 径流 量 的 影响 为 83.04% ,降水 变 
化 对 皇甫 川 流 域 径流 量 影 响 为 16.96%, 人 类 活动 在 
皇甫 川 流 域 径流 量 的 变化 中 起 了 主要 作用 ( 表 3)。 

以 1954 一 1986 年 为 基准 期 ,1987 一 2018 年 人 类 
活动 对 皇甫 川 流 域 输 沙 量 的 影响 为 75.24% ,降水 变 
化 对 皇甫 川 流 域 输 沙 量 影响 为 24.76% ,人 类 活动 在 
皇甫 川 流 域 径 流量 变化 中 起 了 主要 作用 ( 表 4)。 结 
合 皇 惠 川 流域 1954 一 2018 年 降水 量 呈 现 出 不 显著 
下 降 趋 势 ( 图 2) ,而 皇甫 川 流 域 的 水 沙 变化 特征 却 
呈现 出 显 车 减少 趋势 (图 3) ,表明 流域 内 逐年 增强 
的 人 类 活动 是 促使 皇甫 川 流 域 径流 和 泥 沙 减少 的 
主要 因素 ,而 降水 变化 引起 的 降水 减少 是 次 要 


90 
2o; 0 (b) 
£ y=0.004x+4.5285 EN 
2 60 R=0.9971 M. 三 y=0.0013x+0.5418 
ES 1986 年 1 R?-0.996 1986 年 
4 45 EN 
Ts y=0.0004x+15.075 | & y=0.0015x+42.627 
RE R’=0.9098 E R1-0.9178 
ak 15 Ed 
0 ——————— 
5000 10000 15000 20000 25000 30000 5000 10000 15000 20000 25000 30000 
累积 降水 量 /mm 累积 降水 量 /mm 
图 6 降水 量 -径流 量 与 降水 量 - 输 沙 量 双 累 积 曲线 
Fig. 6 Precipitation-runoff and precipitation-sediment transport double cumulative curve 
表 3 皇甫 川 流 域 径流 量变 化 原因 分 析 
Tab.3 Causes analysis of runoff variation in Huangfuchuan basin 
nA sac m Y pe A sa E 实测 年 降水 量 降水 影响 人 类 活动 影响 
时 期 实测 径流 量 计算 径流 量 ep mu) T cU XE — 
7 /108 m? /108 n 减少 里 分 比 减少 量 百分比 减少 量 百分比 
/10* m° /% /10° m° /% /10 m° /% 
1954 一 1986 年 58.13 60.20 - - - - - - 
1987—2018 年 20.12 51.68 38.01 65.38 6.45 16.9696 31.56 83.0496 
表 4 皇甫 川 流 域 输 沙 量变 化 原因 分 析 
Tab.4 Causes analysis of sediment transport variation in Huangfuchuan basin 
实测 年 降水 量 降水 影响 人 类 活动 影响 
实测 输 沙 量 计算 输 沙 量 = l E : — 
时 其 Ah Ho. 减少 量 5 分 比 减少 量 。 ”百分比 减少 量 。 ”百分比 
/10°t /% /10°t /% AO't /% 
1954 一 1986 年 18.33 18.63 - - - - - - 
1987—2018 年 5.TT 15.23 12.56 68.52 3.11 24.7696 9.45 75.2496 
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因素 。 分 析 法 定量 评价 了 降水 和 人 类 活动 对 皇甫 水 文 站 


通过 公式 (2) 和 公式 (3) 分 别 计算 不 同时 期 输 
沙特 征 因子 (S.PY.s) 的 比例 变化 率 ,从 而 确定 各 因 
子 对 皇甫 川 输 沙 量变 化 的 贡献 程度 。 如 表 5 所 示 ， 
1954 一 2018 年 皇甫 川 输 沙 量 总 体 呈 下 降 趋势 , 输 沙 
量 的 比例 变化 率 为 0.49%.a'。 其 中 ,85% 变 化 率 下 


水 沙 变 化 的 贡献 率 ,得 出 降水 量变 化 对 径流 量 的 贡 
献 率 为 25.8% ,降水 量 对 输 沙 量 的 影响 与 之 相似 。 
万 晓 丹 等 采用 线性 趋势 法 对 皇甫 川 流 域 1954 一 
2016 年 径流 量 和 输 沙 量 进行 了 人 研究 ,计算 出 皇甫 川 
流域 多 年 平均 降水 量 为 424 mm, 降 水 量 减少 是 导致 


降 由 产 流 因子 /所 贡献 ,29% 变 化 率 下 降 由 产 沙 因子 
s 所 贡献 ,而 bp 的 变化 趋势 与 输 沙 量 相 反 , 贡 献 了 输 
沙 量 比例 变化 率 的 -4%。 皇 甫 川 流 域 径流 量 和 输 沙 
量 突变 均 是 从 1986 年 开始 ,将 输 沙 量 的 变化 分 为 
1954 一 1986 年 和 1987 一 2018 年 两 个 时 期 进行 分 
析 。 在 1954 一 1986 年 期 间 , 输 沙 量 有 明显 的 上 升 ， 
其 比例 变化 率 为 0.54% .a, 其 中 /因子 贡献 了 比例 
变化 率 的 90% , s 因子 贡献 了 比例 变化 率 的 31% ,而 
尸 因子 贡献 了 比例 变化 率 的 -21%。 在 1986—2018 
年 , 输 沙 量 的 比例 变化 率 为 0.36%.a ,具有 明显 下 
降 趋 势 , 其 中 /因子 贡献 了 输 沙 比例 变化 率 的 91% ,s 
因子 贡献 了 比例 变化 率 的 25% ,P 因 子 贡献 了 比例 
变化 率 的 -16%。 从 3 个 变化 阶段 的 分 析 结 果 可 以 
看 到 ,皇甫 川 流 域 产 流 因 子 / 对 输 沙 量 的 比例 变化 率 
起 着 决定 性 作用 , 产 沙 因 子 ; 对 输 沙 量 比例 变化 率 的 
贡献 相对 较 小 ,而 降水 因子 P 其 至 起 着 相反 作用 。 


表 5 不 同时 期 输 沙特 征 因子 的 比例 
Tab.5 Proportion of sediment transport characteristic 
factors in different periods 


WT 1954—20184E —1954—19864E 1986—2018 4E 
PI% «a^') -0.07 -0.09 -0.06 
Ji (06 a) 0.42 0.49 0.33 
s/(% .am 0.14 0.17 0.09 
总 计 0.49 0.54 0.36 


3 讨论 


降雨 是 径流 补给 的 主要 来 源 ,直接 影响 着 径流 
大 小 ,而 径流 的 变化 会 影响 到 输 沙 变 化 ,所 以 降水 
成 为 影响 径流 和 输 沙 的 决定 因素 之 一 。 和 气象 资料 
研究 表明 ,皇甫 川 流 域 降水 具有 减少 趋势 "5 ,气候 
总 体 呈 暧 干 化 方向 发 展 ,气候 变化 引起 的 降水 减少 
对 皇甫 川 流 域 径流 量 和 输 沙 量 产 生 重要 影响 "。 
1954 一 2018 年 皇甫 川 流 域 降水 呈现 下 降 趋 势 ( 图 
2) ,水 沙 与 降水 有 相同 的 变化 趋势 (图 3) ,表明 降水 
是 影响 流域 水 沙 的 重要 因素 。 赵 广 举 等 ”通过 水 文 


皇甫 川 流 域 水 沙 量 相应 减少 的 重要 因素 。 胡 智 丹 
等 "运用 水 文敏 感性 分 析 法 探究 了 皇甫 川 流 域 径 
流 锐 减 的 原因 ,得 出 皇甫 川 流 域 天 然 时 期 (1985 一 
1998 年) 和 径流 变化 时 期 (1999 一 2006 年 ) 相 比 , 降 
水 减少 对 径流 量变 化 的 贡献 率 为 46%。 以 上 相关 
人 研究 结论 都 与 本 人 研究 结论 一 致 。 但 皇甫 川 流 域 降 
水 变化 呈现 出 不 显著 减 小 趋势 (P>0.1) ,而 水 沙 变化 
却 呈 现 出 显著 趋势 (P<0.05) ,水 沙 下 降 趋 热 相 比 于 
降水 更 加 显著 ,降水 影响 皇甫 川 流 域 产 水 产 沙 能 
随时 间 的 变化 在 减弱 ,并 非 是 影响 水 沙 变化 的 唯一 
素 , 还 有 其 他 重要 因素 对 水 沙 变 化 产生 影响 。 

人 类 活动 是 指 城镇 化 建设 .水土 保持 过 程 E 
产生 活用 水 ,水利 工程 修建 等 各 种 直接 对 流域 径流 
量 和 输 沙 量 或 改变 下 垫 面 条 件 产 生 影响 的 活动 中 。 
皇甫 川 流 域 水 沙 变 化 受 流 域内 一 些 人 类 活动 (退耕 
还 林 还 草 . 淤 地 坝 建 设 .水 保 过 程 等 ) 的 影响 ,使 得 
流域 下 垫 面 发 生 了 明显 改变 ,对 水 沙 扰动 增加 , 减 
少 了 河流 的 径流 量 和 输 沙 量 ,尤其 是 2000 年 皇甫 咱 
流域 开始 了 大 规模 的 退耕 还 林 草 工程 ,造成 流域 内 
大 量 的 耕地 变 为 草地 、 林 地 ,在 一 定 意 义 上 起 到 水 
土 保持 防风 固沙 的 作用 ,其 次 是 一 些 免 耕 . 少 耕 等 
保护 性 耕作 技术 大 力 推广 ,直接 减少 了 水 土 流失 的 
物质 来 源 , 作 物种 植 结构 ,种植 密度 优化 ,使 得 地 表 
裸 土 面积 减少 ,有 效 拦截 降雨 ,减少 地 面 侵蚀 ,改善 
了 土壤 结构 ,使 得 土壤 空 际 增加 ,透水 性 加 强 , 减 少 
了 地 面 径流 和 输 沙 的 效果 所 。 慕 星 等 9 采用 水 文 
分 析 法 定量 评价 了 降水 和 人 类 活动 对 皇甫 水 文 站 
水 沙 变化 的 贡献 率 ,发 现在 20 世 纪 80 年 代 , 人 类 活 
动 对 皇甫 川 流 域 径流 变化 的 贡献 率 为 77.5% ,90 年 
代 人 类 活动 对 皇甫 川 流 域 年 径流 变化 的 贡献 率 为 
56.4%。 刘 微 等 采用 CPA 法 对 皇甫 川 径 流 进行 突 
变 分 析 , 根 据 皇 甫 川 流 域 在 20 世 纪 90 年 代 城 镇 化 
的 发 展 ,城镇 用 地 面积 增 大 ,植被 种 类 改变 ,生态 平 
衡 遭 到 破坏 ,造成 径流 量 随 之 减 小 ,得 出 人 类 活动 
是 2000 年 后 皇甫 川 流 域 径流 量 减少 的 主导 因素 。 
赵 广 举 等 通过 对 皇甫 川 日 径流 量 与 日 输 沙 率 的 变 
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化 特征 研究 发 现 , 在 同期 降水 量 并 没有 显著 减少 的 [2] "iff, FRR, ER. 生态 安全 条 件 下 土地 利用 格局 优化 一 一 以 


情况 下 ,皇甫 川 输 沙 减少 的 主要 原因 是 由 于 河道 内 
存在 的 大 量 洲 地 坝 有 效 地 蓄 水 拦 沙 ,在 1976 年 以 
前 ,皇甫 川 流 域 共 建 成 淤 地 坝 392 座 ,而 到 2015 年 ， 
汶 地 坝 的 数量 达到 了 886 座 ,控制 了 整个 流域 70% 
的 径流 ,河道 径流 被 拦蓄 后 ,上 游 来 水 失去 了 冲刷 
下 游 河道 的 动力 ,从 而 达到 减 刨 的 效果 ,得 出 人 类 
活动 是 导致 皇甫 川 流 域 水 沙 急剧 减少 的 主要 原 
因 。 以 上 相关 研究 结论 与 本 研究 结论 一 致 ,可 以 得 
出 气候 变化 引起 降水 减少 与 人 类 活动 是 导致 皇甫 
川 流 域 径流 量 和 输 沙 量 减少 的 影响 因素 ,其 中 流域 
内 逐年 增强 的 人 类 活动 是 促使 皇甫 川 流 域 径流 和 
输 沙 减 少 的 主要 因素 ,而 气候 变化 引起 的 降水 减少 
是 次 要 因素 。 


4 结论 

本 文 利 用 皇甫 川 流 域 1954 一 2018 年 水 沙 长 时 
间 序 列 监测 数据 ,运用 线性 趋势 法 、Mann-Kendall 检 
验 、 累 积 距 平 法 等 方法 分 析 了 皇甫 川 流 域 水 沙 周 期 
演变 特征 ,应 用 水 文 分 析 法 定量 评估 了 降水 变化 和 
人 类 活动 对 皇甫 川 流 域 径流 量 和 输 沙 量变 化 的 页 
献 率 ,采用 输 沙 归 因 诊断 分 析 判 断 了 降水 、 产 流 能 
力 和 产 沙 能 力 在 皇甫 川 输 沙 量 比例 变化 率 上 的 贡 
献 ,主要 得 出 以 下 结论 : 

(1) 皇甫 川 流 域 降水 量 呈 不 显著 减 小 趋势 ,而 
径流 量 和 输 沙 量 却 呈 显著 减 小 趋势 , 且 径 流量 减少 
趋势 大 于 输 沙 量 的 减少 趋势 。 

(2) 降水 变化 和 人 类 活动 是 造成 皇甫 川 流 域 水 
沙 减 少 的 影响 因素 ,其 中 流域 内 逐年 增强 的 人 类 活 
动 是 促使 皇甫 川 流 域 径流 和 输 沙 减少 的 主要 因素 ， 
而 气候 变化 引起 的 降水 减少 是 次 要 因素 。 

(3) 不 同 阶段 输 沙 归 因 诊断 分 析 表 明 , 产 流 
子 f 对 皇甫 川 输 沙 量 的 比例 变化 率 起 着 决定 性 作 
用 , 产 沙 因子 * 对 输 沙 量 比 例 变 化 率 的 贡献 相对 较 
小 ,而 降水 因子 bp 起 着 相反 作用 。 
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Characteristics of water and sediment variation in the Huangfuchuan basin and 


its influencing factors 
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Abstract: Based on monitoring data of precipitation, runoff, and sediment discharge in the Huangfuchuan basin 
during 1954-2018, the characteristics of water and sediment cycle evolution in Huangfuchuan basin were studied 
using linear trend methods, Mann-Kendall tests, and cumulative anomaly methods. The contribution rates of 
precipitation change and human activities to runoff and sediment transport change in the Huangfuchuan basin 
were quantitatively evaluated using hydrological analysis methods. The results show that precipitation does not 
significantly decrease with time, whereas runoff and sediment transport significantly decrease with time, with a 
greater decreasing trend in runoff than sediment transport. The abrupt change of runoff and sediment transport 
from abundant to dry occurred in 1986. Using 1954-1986 as a reference period, the impacts of human activity on 
runoff and sediment transport in the Huangfuchuan basin from 1987 to 2018 were 83.04% and 75.2495, 
respectively. The impact of precipitation change on runoff and sediment transport in the Huangfuchuan basin was 
16.96% and 24.76% , respectively. Further, reduced runoff capacity contributed more than 85% to changes in 
sediment transport. Precipitation change and human activities are both important factors influencing reduced 
runoff and sediment in the Huangfuchuan basin. Increasing human activities in the basin are the main factors 
reducing runoff and sediment transport in the Huangfuchuan basin, with reduced precipitation caused by climate 
change as a secondary factor. 
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